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Abstract 

In an anisotropic universe, observers not positioned at a point of special symmetry should observe 
cosmic parallax - the relative angular motion of test galaxies over cosmic time. It was recently 
argued that the non-observance of this effect in upcoming precision astrometry missions such 
as Gaia may be used to place strong bounds on the position of off-center observers in a void- 
model universe described by the Lemaitre-Tolman-Bondi metric. We consider the analogous effect 
in anisotropic cosmological models described by an axisymmetric homogeneous Bianchi type I 
metric and discuss whether any observation of cosmic parallax would distinguish between different 
anisotropic evolutions. 
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I. INTRODUCTION 



Homogeneity and isotropy are the cornerstones of the standard cosmological model, pro- 
viding not only a tremendous simplification of General Relativity, but remarkable agreement 
with all observations. Homogeneity is supported by the observed galaxy distribution from 
large-scale structure surveys, while isotropy is supported by, in particular, the remark- 
able uniformity of the average temperature of the cosmic microwave background radiation 
(CMB). Nevertheless, the paradigm-changing observation of cosmic acceleration, now more 
than a decade old, has forced cosmologists to re-examine even their most cherished assump- 
tions, including the correctness of General Relativity, a vanishing cosmological constant, 
and, more recently, homogeneity and isotropy. 

The primary evidence for the accelerating universe comes from the unexpected dimming 



of type la supernovae ], as measured through their light curves. The 

connection to cosmic acceleration requires the assumptions of homogeneity and isotropy, 
and thus this raises the possibility that abandoning one of these principles may allow for 
the appearance of accelerated expansion without actual acceleration itself. 

Of course, the usual cosmological Friedmann- Robert son- Walker (FRW) metric is so sim- 
ple by virtue of its underlying symmetries. Abandoning these leads to a correspondingly 
more complicated metric. It is convenient therefore to begin by studying toy models. One 
class of these that has shown some promise in this direction are the Lemaitre- Tolman-Bondi 
(LTB) metrics, in which we are assumed to live inside a spherically symmetric underdense 
region of spacetime (or "void" ) embedded in an otherwise spatially flat and homogeneous 

I, [isl, Q, [15 1 . Such a spacetime is manifestly inho- 



Einstein-de Sitter universe 



nu je. 
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mogeneous, due to the void, and on its own violates any strong version of the Copernican 
principle, since we must live inside this void in order to account for the observed supernovae 
dimming. Nevertheless, it has been shown that these models can provide a satisfactory fit to 
the luminosity distance-redshift relation of type la supernovae and the position of the first 
peak in the CMB QQ. Thus, LTB models have been suggested as a possible solution 
to the problem of cosmic acceleration, obviating the need for quintessence fields, modifi- 



cations of gravity, or a cosmo^ 



constraining them 



HQ, 
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ogical constant, and considerable effort has been devoted to 



20|. 



Beyond the usual cosmological tests of homogeneity and isotropy, it has recently been 
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suggested that these models could be further constrained by real time cosmological measure- 
ments - in particular, precision measurements of the evolution of the angular positions of 



distant sources. Following the authors of [2l| , we refer to this effect as cosmic parallax. The 
expansion of an FRW universe is isotropic for all observers, and so cosmic parallax would not 
be observed. Of course although our universe is very close to an FRW universe, it is not ex- 
actly so-for example, on large scales bound structures may acquire small peculiar velocities, 
giving rise to a slight deviation from observed isotropic expansion. However, to any observer 
living off-center inside the void of an LTB universe, cosmic evolution itself is anisotropic and 
is an additional source of cosmic parallax. For sufficiently off-center observers the cosmic 
parallax due to anisotropic expansion would dominate over the contribution from peculiar 
velocities. Cosmic parallax could therefore provide an interesting test of void models. Up- 



coming sky surveys such as Gaia 22j may be able to initiate a measurement of this effect, 
requiring only that a similar survey be completed 10 years later in order to complete the 
measurement. The absence of cosmic parallax beyond what is expected from peculiar ve- 
locities would put an upper bound on how far our galaxy could be from the center of the 



void in otherwise allowed LTB models. For example, the authors of [2l| argue that Gaia 
may map sufficiently many quasars, with enough accuracy, so that two such surveys spaced 
10 years apart could detect the additional LTB contribution to cosmic parallax if the Milky 
Way is more than 10 megaparsecs from the center of a 1 Gpc void. If after a decade no 
additional contribution were found, that would constrain the Milky Way to lie unnaturally 
close to the center of such a void. 

On the other hand, detection of a contribution to the cosmic parallax not arising from 
peculiar velocities would indicate that the expansion of the universe is anisotropic from our 
vantage point. In LTB models this would be due to living away from the center of the void, 
and the strength of the additional contribution would be related to this distance. But cosmic 
parallax not attributable to peculiar velocities is a generic feature of any cosmological model 



with anisotropic expansion-an observation also made in 2l|]. For example, spacetimes of the 
Bianchi type would exhibit an additional contribution to cosmic parallax around every point. 
These homogeneous and anisotropi c sp acetimes have recently been invoked during inflation 



to explain anomalies in the CMB 
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25P and during late-time cosmology to describe 



^ Such theories may not, however, be without problems 26[. 
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the expansion driven by anisotropic dark energy |27l . |28j. The fact that motivations exist 
for studying both LTB and Bianchi spacetimes, and that both may exhibit contributions to 
cosmic parallax that are not attributable to peculiar velocities, raises the question of how we 
might interpret any additional cosmic parallax signal. Although an observed signal would 
provide evidence for deviations from an FRW universe, such deviations could in principle be 
due to deviations from spatial homogeneity (as in LTB models) or deviations from isotropy 
(as in Bianchi models). It is therefore of interest to consider how cosmic parallax in Bianchi 
models differs from LTB models. 



The structure of the paper is as follows. In section [TTl we briefly review the results of 21|. 
In section IIIII we describe the kinds of anisotropic models that we will focus on; namely a 
subclass of Bianchi Type I models. Since we know from observations of the CMB that the 
universe is very nearly isotropic at the time of last scattering, we describe how we restrict 
ourselves to Bianchi Type I metrics which pass existing cosmological tests. In section HVl we 
derive the geodesic equations for these spacetimes and discuss our procedures for numerically 
integrating them. We then present the cosmic parallax signal we find for these models and 
discuss how it compares to the cosmic parallax in void models. 

II. COSMIC PARALLAX IN AN LTB VOID 

We begin by briefly reviewing how cosmic parallax manifests itself in LTB models, as 
discussed in [2l|]. Here and throughout the rest of the paper, unless otherwise indicated, by 



cosmic parallax we mean cosmic paral 



The LTB metric is given by 29|, |30|, l31 1 



ax due to anisotropic expansion about an observer. 



ds^ = -de + \^^^^2l)Ldr^ + R\t, r)dn'' , (1) 
1 -|- p(t) 

where R(t,r) is a position-dependent scale factor, /3(r) is related to the curvature of the 
spatial slices, and ()' = d/dr. The Einstein equations relate R(t, r) to [3{r) and an additional 
arbitrary function of integration a{r). Specifying a(r), /3(r), and an initial condition for 
-R(t, r) completely determines the spacetime. In models with an underdense region, or 
"void" , surrounded by an overdense region, a(r) and (3{r) roughly correspond to the width of 
the void and the gradient of the boundary between the inner and outer regions, respectively, 
and will be specified below. 
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The LTB metric describes a region of spacetime that is isotropic about the origin but 
inhomogeneous with respect to the radial direction. Therefore distant galaxies appear to 
be receding at the same rate in all directions for observers located at the origin. On the 



other hand, observers located away from the center o: 



anisotropic recession. One way to observe this effect 2l| is to measure how the angle between 



the void could in principle observe 



the positions of two distant sources evolves over time. This difference is referred to as the 
cosmic parallax. To study this, one considers null geodesies in an LTB universe, obeying 



— + r'' 



d\^ 



""c/A d\ 



. 



(2) 



Here T^pcr are the Christoffel symbols, A is an affine parameter along null geodesies, and 
the four-velocities = ^ satisfy u^u^ = 0. The goal will be to solve these equations to 
determine the null geodesies along which light travels from various sources to an observer, 
and to do this for two different observation times. 



Following 2l|], we work in spherical-polar coordinates {t,r,6,(j)) for which the origin 
coincides with the center of the void (labelled "O" in figured]). Without loss of generality 
we choose the observer to lie on the polar axis at a coordinate distance Tq along it. Spherical 
symmetry about the observer is now broken but the remaining cylindrical symmetry applied 
to (E]) allows us to neglect the polar angle dependence. The system then reduces to three 
second-order geodesic equations, or equivalently six first-order equations. Applying the null 
geodesic condition further reduces the system to five independent first-order equations for 
t, r, 9, p = dr/dX and the redshift z as [l7|, l2l| 



dt 

dX 

dr 

dX 
d9 

dX 

dz 

dX 



p 
J 



[l + z) 



P ^ 



R'R' n R 



dp 
dX 



^ = 2R'P, 



J 

RR' 



1 + 3 , 



R^' 



■r 



R" 



2 + 2/5 R 



(3) 

(4) 
(5) 

(6) 
(7) 



where J = R^d6/dX = Jo, is constant along the geodesic. 
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FIG. 1: In the LTB model the observer is located at a distance rg from the center of the void 
along the axis of symmetry. Each point along the geodesic is described by the spatial coordinates 
(^(A),r(A)) or equivalently by a;*(A) 



To completely specify the system we require five initial conditions, which we provide 
for convenience at the initial observation event, and denote with a subscript "0". Since 
one would like to specify initial conditions in terms of physically measurable quantities, we 
consider the angle between the polar axis and the line of sight along an incoming photon 
trajectory arriving at the observer (see figure [1]). This angle coincides with the coordinate 
angle Q when tq = 0, but in general it is given by 17 1 

i?'(t,r)p 



cos^o 



A=0 



This expression can be used to express Jq and po in terms of to, tq, and via 

Jo 
Po 



(8) 



R{to,ro) sin^o 
cos^o 



Vl+/5(ro) . 



(9) 



Therefore, the system is completely determined by specifying to, tq, ^o, zo, and ^o- 

Clearly, our coordinate choice means that = 0, and our conventional definition of 

I ? 

redshift yields Zo = 0. Following ^] we choose ro = 15Mpc, which is the largest value 



consistent with the CMB dipole [17[ . What remains is to specify a direction on the sky and 
the time of observation to uniquely determine a geodesic that terminates at the spacetime 
event of observation. This picks out an initial geodesic along which light from a distant 
source travels to reach the observer. 

Given the redshift and line of sight angle of a source at an initial time, we can deter- 
mine the trajectory of light from that source at a later time in the following way. As just 
mentioned, the line of sight angle .^o picks out an initial geodesic, terminating at the initial 
observation event. The observed redshift determines how far backwards in time along this 
initial geodesic the source lies. We extract the comoving coordinates of the source by nu- 
merically integrating backwards along this initial geodesic. Once the comoving coordinates 
of the observer and source are determined, and the interval of time between observations 
is specified, we solve a boundary-value problem to determine the final geodesic. We then 
extract the line of sight angle of this final geodesic ^f. 

After repeating this procedure for two sources we have four angles: the input angles ^ao 
and 1^60 and the output angles ^af and ^bf- We may then compute the difference 

^7 = 7/ - 70 = i^af - 6/) - i^aO - 6o) , (10) 

which is the main quantity of interest, hereafter referred to as the cosmic parallax. 



In 2l|] this quantity is calculated for LTB models characterized by the functions 

m = (i70-)%^^(l-tanh^) (11) 



a(r) = (/foOUT^%-3-r/?(r) , (12) 

which characterize the smooth transition from the inner underdense region to the outer, 
higher density region. The quantities Aa, r^o, Ar, and Hq^"^ are free parameters of the 
model that can be tuned to fit CMB and supernovae measurements. Following 2j| and 15[ , 



we choose a model which is in good agreement with SNIa observations and the location of 
the first peak of the CMB, namely Aa = 0.9, r™ = 1.3 Gpc, Ar = 0.4 r™, and H^^^ = 
51 kms~^ Mpc~^. Although this Hubble parameter outside the void seems to be in conflict 



with the measured Hubble Key project value of 72 ± 8kms~^ Mpc~^ 32], one should keep 
in mind that the measurements made to determine this value are all made at less than a 
few hundred megapasec {z ^ 0.01), whereas the size of the void in this model is about 
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1 Gpc {z ^ 0.02). The Hubble parameter inside the void in this model turns out to be 
65kms~^ Mpc~^, which is consistent with the Hubble Key project value. 

The authors of [21 1 numerically compute the cosmic parallax for two different models. 
They find an angle-dependent signal with maximal value roughly 0.18 yuas for sources at 
redshift of 1. It is claimed that this signal is within reach of the forthcoming Gaia mission, 
provided Gaia produces a sky map of the expected number of quasars (roughly 500, 000) 
with final positional error less than 30 /ias, and that such measurements are repeated after 
10 years. 

III. BIANCHI TYPE I MODELS 

As we have discussed, anisotropic expansion of the universe around a given observer 
contributes to cosmic parallax. In the case of LTB models, this may allow us to constrain 
the distance of our galaxy from the center of the void. The further we are from the center 
of the void, the more anisotropic the universe would look to us. Of course, LTB spacetimes 
are not the only ones that can give rise to anisotropic expansion around a point. This raises 
the question of how we might interpret any observation of cosmic parallax that cannot be 
accounted for by peculiar velocities alone, since such an observation would not itself be 
evidence that we live in an LTB universe. To understand this therefore, we study cosmic 



paral 
in 



211. 



ax in alternative anisotropic cosmologies and compare our results with those obtained 



For simplicity we consider anisotropic spacetimes that are spatially homogeneous |33l. 134 1 . 
In particular, we focus on a subclass of Bianchi Type I spacetimes, the metric for which may 
be written as 

ds"^ = -dt^ + aj{t)dx^ + al{t)dy^ + al{t)dz^ . (13) 

In general, ai(t), a2{t)i and a^{t) are independent scale factors, describing how the three 
spatial directions scale with time, which reduces to the standard FRW case when ai(t) = 
02 (t) = 03 (t). We specialize to the case when only one of the scale factors differs from the 
others, say ai{t) = a2{t) 7^ 03 (t), in which case the expansion is axisymmetric. We do so 
because we want to compare our results to observations made by an off-center observer in 
an LTB void. Such an observer will experience axisymmetric cosmic expansion, and so the 
most direct comparison will be to an axisymmetric Bianchi-I universe. There is, however. 
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an important difference in the symmetries around observers in these two spacetimes. In 
the axisymmetric Bianchi-I universe there is an additional plane of symmetry normal to the 
axis of symmetry. The same is not true for an off-center observer in an LTB universe. To 
see this, it is sufficient to consider the extremal case of an observer outside the void, who 
can obviously distinguish the two directions along the polar axis: toward the void and away 
from it. This suggests that cosmic parallax in these two types of anisotropic models will 
differ at least qualitatively, if not in magnitude. 

Setting oi = 02 = ci{t) and 03 = b{t) in (fT3|) . the Einstein equations become 

Hi + 2HaHf, = SttGp (14) 
2Ha + ml = -SuGP, (15) 
Ha + Hi + H, + Hi + HaH, = -8nGP, , (16) 

where an overdot denotes a derivative with respect to t, = Py and are anisotropic 
pressures in the different directions, and we have defined the Hubble parameters Ha = a/a 
and Hh = b/b. The conservation of energy equation in this case is 

p = -2Ha{p + p.) - H,{p + P,) . (17) 

The observational success of FRW cosmology places tight constraints on how anisotropic 
the universe can be. In order to restrict ourselves to solutions that remain close to an FRW 
cosmology, we split each of these exact equations into an FRW part which evolves according 
the FRW equations of motion, and a non-FRW part which we require to remain small, in a 



sense that we will now make clear. A similar approach was used in 35|. We define 



Ha{t) = H{t) + tf{t) (18) 

H,{t) = H{t) + eg{t) (19) 

p(t) = -p{t)+er{t) (20) 

P.(t) = Pyit) = Pit) (21) 

P,(t) = P{t) + es{t) (22) 

where overbars denote the FRW quantities and e is a small perturbative parameter for which 
we will determine an upper bound later. Substituting these definitions into equations (fT^ - 
(fT7|l and collecting powers of e gives the zeroth-order (or background) equations, which are 



just the usual ones of the FRW metric, and the first-order equations 

2H{f + 2g) = SnGr (23) 

6Hf + 2/ = -SttGs (24) 

3H{f + g) + f + g = (25) 

r = -3Hr - p(2/ + g) - Hs (26) 

These constitute four equations in four variables, but only three of these equations are 
independent. To close the system we need additional information, which we obtain by 
assuming an equation of state of the form 

s{t) = ar{t) (27) 

where the parameter a is taken to be constant. Note that this is analogous to P = wp, 
except that it relates the anisotropic component of the pressure to the non-FRW correction 
of the energy density. The value of a will be important in determining whether the non-FRW 
parts of equations (|23l)- (!26l) grow or decay in time. 

Realistic models will be those for which the amount of anisotropy is sufficiently small 
in the past and present. Assuming the anisotropy is set (for example by inflation) to be 
sufficiently small at some early epoch, the question then is whether the anisotropy grows or 
not. In our set-up this corresponds to asking whether the non-FRW parts of equations ( l23|l - 
(l26l) grow or not, and if so, how quickly. We will see that a governs the general behavior of 
the non-FRW quantities, but for a given a, the details will depend on the background (FRW) 
solution. Since we are integrating from the present up to redshifts of order 1, our background 
is well described by the ACDM model. Using this background we can analytically find the 
asymptotic behavior of the non-FRW quantities for different values of the equation state 
parameter a. 

The background energy density and pressure for ACDM are 

P = Pm + pA 
P = -PA, 

where pm is the background energy density of matter and p\ is the effective energy density 
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of the cosmological constant. The background equations become 

3H' = 87rG{pA + Pm) 
2H + 3i?2 = SttGpa 

pm = —'iHpm . 

The solutions can be written simply 



H{t) = Acoth ( -At 



Using ([27D and in 



Pm 

A = 
then gives 



Pa 



sinh -At 
2 



r = — 



(3 + (T)v4coth [-At] +- 



UAttG 



sinh (I At) 



r(t) 



(28) 
(29) 
(30) 

(31) 
(32) 
(33) 

(34) 



which can be integrated to find 



r{t) = C\ cosh ( -At 



sinh -At 
2 



-(3+|f7) 



Using the equation of state (1271) immeditely gives s{t\ Equations 
integrated to find the remaining solutions (for o ^ —9/2) 



A 



fit) = - 

[sinh(|At)]- 

git) = -fit) 



C2 



3ci 



(35) 

and (1231) can then be 

(36) 
(37) 

[sinh (f At)] ^ 

Now we consider the asymptotic behavior of these solutions for t — >■ oo and t — >■ in turn. 
The behavior as t oo is 

'3 



2 [sinh (|At)]^" 



AC2 



r{t) ~ s(t) ~ ci 



sinh -At 
2 



fit) ~ -git) ~ ci 



sinh -At 
2 



2{l+i<T) 

Max[-2,-2(l+^cr)' 



(38) 



(39) 



Here we see that a = —3 is the boundary between growing and decaying solutions. If we 
require all of the non-FRW quantities to decay as t — > oo, then we must restrict ourselves 
to (J > —3. 
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As t — s> 0, the situation is slightly more complicated. The functions f(t) and g(t) both 
contain a piece that depends on the value of a and a piece that does not. The part that 
does not behaves as 

[smh{lAt)Y 
whereas the part that depends on a behaves as 

2[sinh(|At)]'^'+^'^^ 

We require that the expansion history is close to FRW in the far past, which amounts 
to demanding that \f /H\, \g/H\, etc remain < 0{1) as t 0. From the first term, this 
requires that C2 = 0, while the second term requires that a < —3/2. Therefore solutions with 
decaying anisotropy have an equation of state parameter lying in the range —3 < a < —3/2. 
Turning to r{t) and as t ^ we have 



r(t) ~ s(t) ~ Cl 



sinh ( -At 
2 



1 -i^+h) 



(40) 



When a > —3 these solutions diverge as t — > 0, but if a < —3/2 they diverge slower than p 
and P, respectively. So the condition 

-3<a<-^ (41) 

will ensure that all non-FRW quantities remain small in the far past, as required. We will 
restrict ourselves to these models in the rest of the paper. 

What remains is to fix ci and the constant A in the solutions (|55]) - (|57|) . We do this by 
imposing observational constraints on the models of interest. First we impose a condition at 
the surface of last scattering. For an anisotropically expanding universe to be viable it must 
at the very least predict angular variations in the temperature of the CMB no bigger than 
10^^. We can estimate the maximum temperature difference at the time of last scattering 
by 

Tiss 1 — V Tiss f ^1 , (42) 



where is the temperature along the axis of symmetry and T^^ is the temperature in the 
transverse plane; subscripts "0" and "Iss" refer to quantities today and at the last scattering 
surface, respectively. Recall that b{t) is the scale factor along the axis of symmetry and a{t) 
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is the scale factor in the transverse plane. From the definitions of Ha and if;,, we have 

a{t) = aoexp [ {H{t') + ef{t'))dt' (43) 

J to 

h{t) = boexp [ {H{t') + tg{t'))dt' . (44) 

J to 

Using these, and choosing oq = &o = 1) we can re- write the expression for ATq to first order 
in e as 



ATn 



rp I -'ISS 

J-lss I ^ 



b\ss cto 



f{t)dt 



to 



(45) 



where we have used the fact that g{t) 
Inserting the solution for f{t) gives 



-/(t), which is a consequence of requiring C2 = 0. 



ATq 



rtiss 






sinh Q^^^ 


'to 





-2(1+1) 



dt 



(46) 



36] yields one condition on the 



Demanding the left-hand side to be at most 1.3 x 10' 
product eci and A for a given equation of state parameter a. 

In order to break the degeneracy between eci and A, we need one further condition. 



Equation ( 1461) already requires the difference between Ha and H^ to be small - well within 



32|. We find it 



the accepted uncertainty in the measured value of the Hubble parameter 
convenient to choose to set the arithmetic average of the Hubble parameters in the three 
directions equal to the observed value. 



(47) 



where again we have used the fact that when C2 = 0, g{t) = —f{t). Alternative choices, such 
as setting Ha or H^, equal to the measured value of the Hubble parameter, would not change 
the order of magnitude of eci and A and would leave the final result essentially unaltered. 
Inserting the solution for f{t) into equation fj471) gives a second condition on eci and A 

2(1+1) 



eci 



A 



[Hohs — H) 



sinh —At 
2 



(48) 



We can (numerically) solve this equation for A in terms of eci and then substitute it back 
into the first condition (I46p to obtain an upper bound on eci. By taking the maximal allowed 
value for eci, we can then find A using (HHl) . For a = —2 we find A ~ 62 and eci ~ 1.3 x 10^^. 
We always take Ci to be C(l), and so in this case we choose A = 1.3 and e = 10~^. 
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In this way we can fully determine solutions to the non-FRW quantities for a given 
equation of state parameter o. By restricting ourselves to —3 < o < —'i/2 we have chosen 
to focus on models for which these solutions decay in the distant future and which diverge 
slower than the respective background quantities in the far past as one approaches the 
initial sigularity. These models seem to be the most conservative realizations of a Bianchi-I 
cosmology, in the sense that they are the easiest to make consistent with observations, or 
alternatively, the most difficult to rule out. One might try to push the boundaries slightly, 
for example by considering models with anisotropics that approach a constant in the distant 
future rather than vanishing. It may be that such models can be carefully tuned to match 
observations. We do not consider these more general models here, since our main interest 
is not model-building, but rather to explore a general effect (cosmic parallax) arising in a 
Bianchi-I universe. 



IV. GEODESICS AND COSMIC PARALLAX IN AXISYMMETRIC BIANCHI-I 
MODELS 

As in section [Tll in order to analyze cosmic parallax we need to find null geodesies in 
the spacetime that connect an observer and various sources at two different times. Letting 
latin indices run over 1, 2, 3 and x^, x^, x^) = {t, x, y, z), the non-zero Christoffel symbol 
components for the Bianchi-I metric in ([T3l) are 

where no sum on the index i is implied and Hi is defined as above. The four geodesic 
equations are then 



i ^ 

(i^x* dt dx' 



-E^. (49) 



dX^ d\ d\ 



2^^:7T3r ' (50) 



with the null geodesic condition, u^Uf^ = 0, becoming 



dX J ^ V d\ 

As before, = dx^^/dX. We specialize to the axisymmetric case by setting ai = 02 = a{t) 
and as = b{t), and also Hi = H2 = Ha and H3 = Hb. The scale factors, a{t) and b{t), and 
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Hubble parameters, Ha{t) and Hb{t), are fixed after clioosing a and solving the full set of 
equations as in the previous section. 

After fixing the background Bianchi-I spacetime, equations fl49l) and flSOp yield four 
second-order differential equations and one constraint equation in four dependent variables. 
To solve this system we must in principle specify initial conditions for the four dependent 
variables as well as initial velocities (derivatives with respect to A) giving a total of eight 
conditions. However, using the constraint equation the system can be reduced to seven 
independent first-order equations. 

Considering the u^{t) = dx^(t)/dX as functions of time along the geodesic, equations (l50ll 
can be integrated immediately to give 

= ular\t), (52) 
which can then be used in equation (15T|) to give 

As in the LTB case it is again useful to find an expression for the redshift as it will be one of 
our observational inputs. To find this expression, as usual, we consider two photons emitted 
from a source at times t and t + t, respectively, where r is taken to be infinitesmally small. 
The trajectory of the first photon is described by equation fl53l) . while to first order in r, 
the trajectory of the second photon is described by the geodesic equation 

^ ' i 

Here the variation in time corresponds to a time-delay, and so to a change in geodesic, and 
not to a change of the time coordinate along a fixed geodesic. Since the m*(A) are directional 
derivates along a given geodesic, they remain unaffected by time variations (derivatives) in 
obtainig this equation. Using the standard definition of redshift, l + z{Xem) = T(Aob)/T(Aem), 
we find the relation 

d\og{l + 2:(Acm)) ^ 1 C?r(Acm) ^gg^ 

dXem T(Acm) C^Acm 

where Aem and Aob are the values of the affine parameter at the emission and observation 
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event, respectively. Using this in (IMI) gives 



rflog(l + z) 



dX 



J2i {u'aif Hi 



d\og^J^~Jvht~f 



(56) 



dX 



In the last step we have also made use of fl5^ and the chain rule 



d 

dX 



{u\X)ai{X)) 



dt 
dX 



{u'{t)a,{t)Hi{t)) . 



With the initial condition z{Xq) = 0, integrating fISB]) then gives 



(1 + ^) 



(57) 



As in section HIl we would like to express our results not only in terms of redshift but also 
in terms of angles, since these are the actual observables. In general there are four pieces 
of data for each object in the sky, namely the time of observation, two angles with respect 
to an arbitrary coordinate system, and the observed redshift of the source. In the case 
of cylindrical symmetry one angle is enough. In contrast with LTB spacetimes, not much 
is gained in this case by rewriting the geodesic equations in terms of angles and redshift. 
Instead, we numerically integrate the equations in the above coordinates and then express 
the results in terms of angles and redshift. 

Our procedure for computing the cosmic parallax for these models is analogous to that 
in section HIl We work in local Cartesian coordinates (t, Xi, 0:2, Xs) in which the observer is 
located at the origin. As seen above, the system of four second-order differential equations 
plus a constraint reduces to seven independent first-order equations, some of which can 
be integrated immediately by hand. We are then left with the problem of fixing initial 
conditions for our complete set of equations 



u%X) 




(58) 



(59) 



dx'^ 

Ix 



(60) 
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where the subscripts o refer to quantities at the initial observation event, corresponding to 
A = 0. To close the system requires seven initial conditions that unfortunately cannot be all 
specified at the observation event. By construction we have x*o = 0, leaving four remaining 
conditions, three of which are obtained by specifying and two initial spatial velocities in 
terms of measurable quantities (angles), while the last one is given by the observed redshift. 
We integrate backwards along the initial geodesic until reaching the desired redshift and then 
find the comoving coordinates of the source. To find the final geodesic that connects the same 
source with the observer at a later time to + At, we solve the corresponding boundary- value 
problem (namely, to find solutions of the null geodesic equations with two fixed endpoints). 
We then find the velocities along the final geodesies at the time of observation. 

This procedure is repeated for two sources, yielding four sets of velocities: one set for 
each initial geodesic and one set for each final geodesic. The velocities and f * along two 
geodesies at the same observing time (see figure [2]) are related to the angle 7 between them 
by 

cos (7) = ^-"-"'"^ . (61) 



.k\2 



We then calculate the cosmic parallax A7, as defined in (fTOl) . by taking the difference of the 




XY 



FIG. 2: Initial angle 70 defined by the velocity vectors at observing time for two sources 



angle between the two sources at the two different times. Finally, we plot the cosmic parallax 
as a function of the polar coordinate 6 for one of the sources. By convention we choose the 
second, or trailing, source as our reference, as shown in figure [31 We find the parallax for 
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z 




XY 



FIG. 3: Location of sources as seen by an observer at the measuring event. The vertical axis points 
along the axis of symmetry. We are considering sources at equal redshifts of z = 1 and in a plane 
defined by a fixed value of the polar angle (f). 

two hypothetical sources with the same redshift, initially separated by a given angle on the 
sky, and for a given At and e. In figure H] we plot the parallax for two hypothetical sources 
at 2; = 1, with an initial separation of 90 degrees on the sky, for e = 10^^ and various values 
of At. For At = 500 years and e = 10~^ we find this maximal value to be of the order 
of 6 X 10~^^ radians. Of course, our true goal is to find the maximal signal for reasonable 
time scales, say At ~ 10 years, but unfortunately numerical noise dominates over the signal 
for At of that magnitude. This problem is shown in figure H] where the results from our 
code for the cosmic parallax are plotted for different values of At. One can see that by 
At = 20 years the signal-to-noise ratio becomes very poor, making it difficult to extract 
trustworthy predictions. Therefore, although we can directly compute the signal for this 
model at At = 10 years, we prefer to calculate the cosmic parallax for several values of At 
between 5 and 500 years and interpolate the value at 10 years. Since our primary goal is 
to find an order of magnitude estimate for the effect, this approach should be acceptable. 
Figure [5] shows the values of the cosmic parallax between two sources for decreasing values 
of At (keeping all other values fixed). For each value of 9, we use a linear fit passing through 
the origin (because cosmic parallax must vanish for At = 0) to find the interpolated value 
at 10 years. We could just as well perform an extrapolation by omitting data below some 
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FIG. 4: A sequence of cosmic parallax signals for different values of At. Top row, left to right: 
At = 20yrs, At = 80yrs. Bottom row, left to right: At = 120yrs, At = 200yrs. The signal-to-noise 
ratio becomes smaller as At becomes smaller. 



cut-off, say At = 50 years. This does not have any noticeable effects on our results. 

After repeating the interpolation for all of our data points, the result is plotted in figure El 
In the same figure we also fit a sine function through the parallax data obtained directly from 
our code for At = 10 years. This fit is quite close to the interpolated values, suggesting that 
the numerical noise is randomly distributed about zero and strengthening our confidence on 
the correctness of the linear interpolation procedure. 

For the largest allowed values of e consistent with the observed anisotropy of the CMB 
(i.e., e ~ 10~^) and for time delays on the order of 10 years, we find that the maximal cosmic 
parallax is on the order of 10~^^ radians, or equivalently 10"'' microarcseconds. This is three 



orders of magnitude smaller than the maximal cosmic parallax seen by 



2l| for LTB models. 



It is also three orders of magnitude smaller than the expected level of cosmic parallax from 



peculiar velocities alone in a ACDM universe 



371] . In other words, the contribution to cosmic 



parallax in this model due to anisotropic expansion is sub-dominant to the contribution from 
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Ayifm) 



371 

4 



FIG. 5: Cosmic parallax as a function of At. Here e = 10 ^ and sources are at z = 1, separation 
angle of 90 degrees. 



peculiar velocities ^. 

The qualitative behavior of the cosmic parallax in this model is also quite different from 



that of LTB models as put forth in [2l|]. Both LTB and Bianchi models show a sinusoidal 
(or at least quasi-sinusoidal) cosmic parallax. However, whereas LTB models exhibit a 2it- 
periodic behavior, here we see that Bianchi models exhibit a vr-periodic behavior. This is 
to be expected due to the symmetries of the two types of spacetimes, as alluded to earlier 
in the paper. The LTB spacetime is axisymmetric about an off-center observer but is not 
plane-symmetric about the plane normal to the symmetry axis. If we align the z-axis along 
the symmetry axis, then in spherical-polar coordinates this amounts to saying that the 
spacetime is invariant under changes in (the azimuthal angle) but has no symmetries 
under (nontrivial) changes in 6 (the polar angle). In other words, one would expect cosmic 
expansion to be 27r-periodic in 6, which is just what we see for the cosmic parallax in these 
models. The Bianchi spacetimes, on the other hand, are both axisymmetric and plane- 
symmetric about the plane normal to the axis of symmetry. So one would expect cosmic 



^ Here we assume that the Bianchi-I models we consider have roughly the same peculiar velocity-redshift 
relation as an FRW universe 
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At= lOyr 

Ay(/ias) 




e(rad) 



FIG. 6: Estimate of the cosmic parallax for At = 10 years, e = 10^ , and for sources at z = 1 with 
a separation angle of 90 degrees, plotted as a function of the angle between the trailing source and 
the symmetry axis (i.e., the angle between the leading source and the symmetry axis is advanced 
by 90 degrees). Each of the plotted points is an interpolated value using 100 trials at values of At 
spread from 5 — 500 years. The solid line is a fit of a sine function through the raw data for the 
cosmic parallax calculated for At = 10 years. 

expansion to be vr-periodic in 9, which is what we see for the cosmic parallax in these models. 

Although here we have only considered a particular model, preliminary investigations into 
other models suggest that these results are robust. For example, one might consider models 
with other constant values for the equation of state parameter o", a time- varying rather than 
constant cr, or an equation of state that takes a different form than equation (1271) . None 
of these modifications seem to affect the order of magnitude of the cosmic parallax (which, 
due to the symmetry of the metric, is the only free parameter). What this suggests is that 
the contribution to the cosmic parallax from viable Bianchi-I models is much smaller than 
the contribution from viable LTB-void models. If the observed cosmic parallax deviates 
from what is expected in an FRW universe, it is unlikely that this is due to our living in a 
Bianchi-I spacetime. 
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V. CONCLUSIONS 



In order for the standard FRW cosmology to agree with observations, the expansion of 
the universe must be accelerating. So far, all suggested mechanisms to drive such accel- 
eration involve new physics; either the existence of exotic new components of the cosmic 
energy budget, modifications to Einstein's theory of gravity, or a cosmological constant. 
Alternatively, it has been suggested that the phenomenon of cosmic acceleration is due to 
interpreting results in the FRW model, when in fact the correct underlying cosmic geometry 
could be that of a void model, such as that idealized by an LTB spacetime. These models are 
relatively simple and make a host of predictions that can be used to either test or constrain 
them. One prediction that may soon be testable is cosmic parallax. For observers located 10 
Mpc from the center of a 1 Gpc LTB void, this effect would have a magnitude of the order 
of 10~^/ias per decade for sources at redshift 1. Since the effect in an FRW universe (due to 
peculiar velocities) is expected to be roughly the same order of magnitude, one might hope 
to subtract the signal due to peculiar velocities from the signal due to anisotropic expansion 
about a point. Measuring this additional contribution to the cosmic parallax would clearly 
indicate a departure from FRW. 

However, what is less clear is how we might interpret detection of an additional con- 
tribution to cosmic parallax. Here we have examined an axisymmetric Bianchi-I universe 
as an alternative explanation for any cosmic parallax component that is not attributable 
to peculiar velocities. We find that for a class of models whose anisotropy is consistent 
with the observed temperature anisotropies of the CMB and Hubble expansion today, the 
maximum amount of cosmic parallax is three orders of magnitue smaller than the maximal 
signal in LTB models. Perhaps more importantly, the maximum effect is also three orders 
of magnitude smaller than the expected level of cosmic parallax from peculiar velocities in 
an FRW universe. Although we have focused our discussion in this paper on a particular 
model, we have found these results to be fairly model- independent. Thus it seems unlikely 
that measurements of cosmic parallax can constrain Bianchi-I models that are not already 
ruled out by the CMB or other cosmological observations. Therefore, while cosmic parallax 
will be nonzero for any anisotropic expansion, the magnitude of the effect suggested in 
appears to be significantly larger, and qualitatively different than in the class of models 
considered here. 



2l| 
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